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Abstract 

Photoproduction of beauty quarks in events with two jets and an electron asso- 
ciated with one of the jets has been studied with the ZEUS detector at HERA 
using an integrated luminosity of 120 pb~^. The fractions of events containing b 
quarks, and also of events containing c quarks, were extracted from a likelihood 
fit using variables sensitive to electron identification as well as to semileptonic 
decays. Total and differential cross sections for beauty and charm production 
were measured and compared with next-to-leading-order QCD calculations and 
Monte Carlo models. 
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1 Introduction 



The production of heavy quarks in ep coUisions at HERA is an important testing ground 
for perturbative Quantum Chromodynamics (pQCD) since the large 6-quark and c-quark 
masses provide a hard scale that allows perturbative calculations. When Q^, the negative 
squared four-momentum exchanged at the electron or positron^ vertex, is small, the re- 
actions ep —>■ ebbX and ep ^ eccX can be considered as a photoproduction process in 
which a quasi-real photon, emitted by the incoming electron interacts with the proton. 

The corresponding leading-order (LO) QCD processes are the direct-photon process, in 
which the quasi-real photon enters directly in the hard interaction, and the resolved- 
photon process, in which the photon acts as a source of partons which take part in the 
hard interaction. For heavy-quark transverse momenta comparable to the quark mass, 
next-to-leading-order (NLO) QCD calculations in which the massive quark is generated 
dynamically [1,2] are expected to provide reliable predictions for the photoproduction 
cross sections. 

Beauty and charm quark production cross sections have been measured using several 
different methods by both the ZEUS [3-18] and the HI [19-30] collaborations. Both the 
deep inelastic scattering (DIS) and photoproduction measurements are reasonably well 
described by NLO QCD predictions. 

Most of the previous measurements of 6-quark production used muons to tag semileptonic 
decays of the B hadrons. The identification of electrons close to jets is more difficult 
than for muons, but the electrons can be identified down to lower momenta. A first 
measurement of 6-quark photoproduction from semileptonic decays to electrons (e~) was 
presented in a previous publication [6], which used the e~^p collision data from the 1996- 
1997 running period corresponding to an integrated luminosity of 38pb~^. This paper 
presents an extension of this measurement exploiting semileptonic decays to electrons as 
well as to positrons for data taken with both e~p and e^p collisions using three times 
the integrated luminosity. The production of electrons from semileptonic decays (csl), 
in events with at least two jets (jj) in photoproduction, ep ebb X ejjeshX', was 
measured in the kinematic range 

< 1 GeV^ and 140 GeV < W^p < 280 GeV, where 
W^p is the centre-of-mass energy of the 7p system. The likelihood method used to extract 
the 6-quark cross sections also allowed the corresponding c-quark cross sections to be 
extracted. This paper provides a complementary study to the measurements using muon 
decays. 



Hereafter unless explicitly stated both electrons and positrons are referred to as electrons. 
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2 Experimental set-up 



This analysis was performed with data taken from 1996 to 2000, when HERA colhded 
electrons or positrons with energy = 27.5 GeV with protons of energy Ep = 820 GeV 
(1996-1997) or 920 GeV(1998-2000). The corresponding integrated luminosities are 38.6± 
0.6 pb~^ at centre-of-mass energy y/s = 300 GeV, and 81.6 ± 1.8 pb~^ at y/s = 318 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [31]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [32], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consisted 
of 72 cylindrical drift chamber layers, organised in nine superlayers covering the polar- 
angle^ region 15° < 6 < 164°. The transverse-momentum resolution for full-length tracks 
is a{pT)/pT = 0.0058pr © 0.0065 © 0.0014/pt, with pt in GeV. The pulse height of the 
sense wires was read out in order to estimate the ionisation energy loss per unit length, 
dE/dx (see Section 3). 

The high-resolution uranium-scintillator calorimeter (CAL) [33] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a{E)/E = 0.18/V^ for electrons and 
a{E)/E = 0.35//E for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — *• ejp, 
where the photon was measured in a lead-scintillator calorimeter [34] placed in the HERA 
tunnel at Z = —107 m. 

3 dE/dx Measurement 

A central tool for this analysis was the dE/dx measurement from the CTD. The pulse 
height of the signals on the sense wires was used to measure the specific ionisation. This 
pulse height was corrected for a number of effects [35]. Such as a factor 1/ sin 6* due to 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj — — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. The azimuthal angle, 0, is measured with respect to the X axis. 
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the projection of the track onto the direction of the signal wire, the space-charge effect 
caused by the overlap of the ionisation clouds in the avalanche, and the dependence of 
the pulse shape on the track topology. An additional correction was needed for hits close 
to the end-plates of the CTD. If a hit followed a previous one on the same wire within 
100 ns, its pulse could be distorted: such hits were rejected. The event topology was used 
to identify additional double hits that could not be resolved; the dE/dx measurement 
was corrected accordingly. 

The dE / dx value of a track was calculated as the truncated mean value of the individual 
measurements, corrected as discussed above, after rejecting the lowest 10% and the highest 
30% of the measurements. Hits where the measured pulse height was in saturation were 
always rejected in forming the mean. Corrections were applied for the finite number 
of hits and whenever more than 30% of the hits were saturated. The corrected dE/dx 
measurement was normalised in units of mip (minimum ionising particles) such that the 
minimum of the dE/dx distribution was 1.0 mip. Electrons are expected to have a mean 
value of about 1.4 mip in the momentum range studied here. 

Different samples of identified particles were used to calibrate and validate the dE/dx 
measurement. The samples used for calibration were: 

• from photon conversions, J/ip decays and DIS electrons; 

• vr^ from decays with 0.4 GeV < p < 1 GeV, where p is the measured track momen- 
tum. 

The samples used for validation were: 

• vr^ from outside the momentum range used for the calibration sample, 
as well as vr^ from A and D* decays; 

• K"^ from 0" and D* decays; 

• p,p from A decays; 

• cosmic fi^. 

Typical sample purities were above 99% for the calibration samples and well above 95% 
for the validation samples [35]. 

After all corrections, the measured dE/dx depended only on the ratio of the particle's 
momentum to its mass, jSj. This is illustrated in Fig. 1. It shows the specific energy loss 
as a function of jS'j, for the different samples of identified particles, e^, fi^,7i^, K^,p,p. 
All particle types are well described using a single physically motivated parametrisation 
of the mean energy loss as a function of f3'y with five free parameters following Allison 
and Cobb (AC) [36]. 
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Given the quality of the description of the mean dE/dx by the AC parametrisation, 
the measurements can be used to determine residuals on dE/dx. As an example, the 
distribution of residuals for a sample of tracks with the number of hits after truncation, 
'T'trunc, equal to 23 is shown in Fig. 2. The dE/dx resolution is typically 11% for tracks that 
pass at least five superlayers. It improves to about 9% for tracks that pass all superlayers. 

4 Monte Carlo simulation 

To evaluate the detector acceptance and to provide the signal and background distri- 
butions, Monte Carlo samples of beauty, charm and light-flavour events generated with 
Pythia 6.2 [37] were used. 

The production of 66-pairs was simulated following the standard Pythia prescription 
with the following subprocesses [38]: 

• direct and resolved photoproduction with a leading-order massive matrix element; 

• b excitation in both the proton and the photon with a leading-order massless matrix 
element. 

The CTEQ4L [39] parton distributions were used for the proton, while GRV-G LO [40] was 
used for the photon. The 6-quark mass parameter was set to 4.75 GeV. The production of 
charm and light quarks was simulated for both direct and non-direct photoproduction with 
leading-order matrix elements in the massless scheme using the same parton distributions 
as for the bb samples. 

The generated events were passed through a full simulation of the ZEUS detector based 
on Geant 3.13 [41]. The ionisation loss in the CTD was treated separately using a 
parametrisation of the measured data distributions based on the calibration sample [38, 
42]. The final Monte Carlo events had to fulfil the same trigger requirements and pass 
the same reconstruction programme as the data. 

5 Data selection 

Events were selected online with a three-level trigger [31,43] which required two jets 
reconstructed in the calorimeter. 

The hadronic system (including the decay electron) was reconstructed from energy-flow 
objects (EFOs) [44] which combine the information from calorimetry and tracking, cor- 
rected for energy loss in inactive material. Each EEC was assigned a reconstructed four- 
momentum g* = [p^^ , Py , p^^ , E^) , assuming the pion mass. Jets were reconstructed from 
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EFOs using the algorithm [45] in the longitudinally invariant mode with the massive 
recombination scheme [46] in which q^^^ = ^ . and the sum runs over all EFOs. The 
transverse energy of the jet was defined as E^r^^ = E^°^ ■ p^^^/p""*, where E^°^, p^"^* and p^^* 
are the energy, momentum and transverse momentum of the jet, respectively. The trans- 
verse energy, -E:^*, is therefore always larger than the transverse momentum, pip", used in 
a previous publication [5]. 

Dijet events were selected as follows: 

• at least two jets with E}^^ > 7(6) GeV for the highest (second highest) energetic jet 
and pseudorapidity of both jets \ri^^^\ < 2.5; 

• the Z coordinate of the reconstructed primary vertex within l^ytxl < 50 cm; 

• 0.2 < yjB < 0.8, where i/jb = {E — Pz)/{2Ee) is the Jacquet-Blondel estimator [47] 
for the inelasticity, y, and E — Pz = Yli ~ P^zi where the sum runs over all EFOs; 

• no scattered-electron candidate found in the calorimeter with energy E'^ > 5 GeV and 
He < 0.9, with He = 1 — (1 — cos6'g), where Q'^ is the polar angle of the outgoing 
electron. 

These cuts suppress background from high-Q^ events and from non-ep interactions, and 
correspond to an effective cut of < 1 GeV^. 

6 Identification of electrons from semileptonic decays 

Electron candidates were selected among the EFOs by requiring tracks fitted to the pri- 
mary vertex and having a transverse momentum, p|n, of at least 0.9 GeV in the pseu- 
dorapidity range \ff\ < 1.5. Only the EFOs consisting of a track matched to a single 
calorimetric cluster were used. To reduce the hadronic background and improve the 
overall description, at least 90% of the EFO energy had to be deposited in the electro- 
magnetic part of the calorimeter. Electron candidates were required to have a track with 
'^trunc > 12 to ensure a reliable dE/dx measurement. An additional preselection cut of 
dE/dx > 1.1 mip was applied to reduce the background. Candidates in the angular re- 
gion corresponding to the gaps between FCAL and BCAL as well as between RCAL and 
BCAL were removed using a cut on the EFO position [48]. 

Electrons from photon conversions were tagged and rejected based on the distance of 
closest approach of a pair of oppositely charged tracks to each other in the plane perpen- 
dicular to the beam axis and on their invariant mass [6] . Untagged conversion background 
and electrons from Dalitz decays were estimated from Monte Carlo studies. 

The electron candidate was required to be associated to a jet using the following proce- 
dure: 
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the jet had to have E^^* > 6 GeV and < 2.5; 



• the distance AR = ^/(r/j^* - r/^)^ + - 0^)^ < 1.5; 

• in case of more than one candidate jet, the jet closest in AR was chosen. 

For the identification of electrons from semileptonic heavy-quark decays, variables for 
particle identification were combined with event-based information characteristic of heavy- 
quark production. For a given hypothesis of particle, i, and source j, the likelihood, Cij, 
is given by 

I 

where 'Pij{di) is the probability to observe particle i from source j with value di of a 
discriminant variable. The particle hypotheses i G {e, ^,n, K,p} and sources, j, for 
electrons from semileptonic beauty, charm decays and background, j G {6, c, Bkg}, were 
considered. For the likelihood ratio test, the test function, Tij was defined as 

~ a a' C ■ 

m,n 

The ctj, a'j denote the prior probabilities taken from Monte Carlo. In the sum, m, n run 
over all particle types and sources defined above. In the following, T is always taken 
to be the likelihood ratio for an electron originating from a semileptonic 6-quark decay: 
T = Te^bi unless otherwise stated. The following five discriminant variables were combined 
in the likelihood test: 

• d-E/dx, the average energy loss per unit length of the track in the CTD; 

• £'emc/-£'CAl, the fraction of the EFO energy taken from the calorimeter information, 
EcAL, which is deposited in the electromagnetic part, E^mc] 

• -EcAL/Ptrack: the EFO energy divided by the track momentum. 

In order to distinguish between electrons from semileptonic 6-quark and c-quark decays 
and other electron candidates, the following additional observables were used: 

• p^', the transverse-momentum component of the electron candidate relative to the 
direction of the associated jet defined as 

rcl _ IPjet X Pel 
P± - 



Wietl 

where pe is the momentum of the electron candidate. The variable p^' can be used 
to discriminate between electrons from semileptonic 6-quark decays and from other 
sources, because its distribution depends on the mass of the decaying particle. It is 
not possible to distinguish charm from light-fiavour decays with this variable; 
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• A0, the difference of azimuthal angles of tlie electron candidate and the missing trans- 
verse momentum vector defined as 

A0=|0(pe)-0(^T)|, 

where is the negative vector sum of the EFO momentum transverse to the beam 
axis, 

and the sum runs over all EFOs. The vector is used as an estimator of the di- 
rection of the neutrino from the semileptonic decay. The variable A0 can be used to 
discriminate semileptonic decays of b quarks and c quarks from other sources. 

The shapes of the charm- and light-quark p^' distributions in the Monte Carlo were 
corrected [5] using a dedicated background sample in the data. The value of the correction 
increased with p^' and was 15% at p'j^^ = 1.5 GeV, where the purity of the b contribution 
is highest. For the A0 distribution a correction was determined in a similar way, but in 
this case the maximal correction was only of the order of 5%. 

In Fig. 3 the distributions of the five input variables used in the likelihood are shown for 
electrons from 6-quark and c-quark decays and for electron candidates from other sources. 
A clear difference in shape between signal and background can be seen. 



7 Signal extraction 

The electron candidates in the Monte Carlo samples were classified as originating from 
beauty, charm or background. The beauty sample also contains the cascade decays b — >■ 
c ^ e, but not 6 — > r — > e and b J/ip ^ e^e~ that have been included in the 
background sample. Test functions (see Section 6) were calculated separately for the 
three samples. The fractions of the three samples in the data, f^^^^^, f^^c^^y fEkg^y were 
obtained from a three- component maximum-likelihood fit [49] to the T distributions. The 
constraint f^^^^^ + ff^f^^ + fEkg^ — ^ imposed in the fit. The fit range of the test 
function was restricted to — 21nT < 10 to remove the region dominated by background 
and where the test function falls rapidly. The for the fit is x^/iidf = 13/12 and the 
6-quark and c-quark measurements have a correlation coefficient of —0.6. The result of 
the fit is shown in Fig. 4 and corresponds to a scaling of the cross section predicted by 
the beauty Monte Carlo by a factor of 1.75 ± 0.16 and the charm Monte Carlo by a factor 
of 1.28 ±0.13. These factors are applied to Figs. 5-8 and denoted as "PYTHIA (scaled)". 
A fit over the whole T range gave consistent cross sections and was used as a cross-check. 

The distributions of the five variables that entered the likelihood are shown in Fig. 5. 
The description of all variables is reasonable. These distributions are dominated by the 
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background contribution. In order to select a beauty-enriched sample, a cut of —2 InT < 1 
was applied. The resulting distributions are shown in Fig. 6. A likelihood for semileptonic 
charm can also be constructed, Tg^c- The distributions of the likelihood for a sample 
satisfying —2 In Tg^c < 1-5 are shown in Fig. 7. Good agreement is observed in both cases. 

To demonstrate the quality of the data description by the Monte Carlo, the distributions of 
E"^* and r]^'^^ of the jet associated with the electron and of the of the electron candidates 
are compared in Fig. 8a)-c). In Fig. 8d)-i) the same distributions are compared for the 
beauty- and charm enriched-samples. Some differences are observed in the jet variables, 
mainly in the region dominated by background. The agreement significantly improves for 
samples enriched in beauty and charm signal. 



8 Cross section determination 

The cross sections have been measured in the kinematic range < 1 GeV^, 0.2 < y < 0.8, 
with at least two jets with E}^^ > 7(6) GeV, \r]^'^^\ < 2.5 and an electron from a semileptonic 
decay with > 0.9 GeV in the range \ri'^\ < 1.5. 

The differential beauty cross section for a variable, v, was determined separately for each 
bin, k, from the relative fractions in the data obtained from the fit and the acceptance 
correction, A^^, calculated using the Monte Carlo, 

dvk Al^ ■ C ■ Avk 

where _/V^^™ is the number of electron candidates found in the data, C is the integrated 
luminosity and Av^ is the bin width. 

In order to determine the acceptance, the jet-finding algorithm was applied to the MC 
events after the detector simulation and at hadron level. The acceptance is defined as 

yyobs 

A = ^ 

where A''°'^^ is the number of electrons from semileptonic decays reconstructed in the 
Monte Carlo satisfying the selection criteria detailed in Sections 5 and 6, and A'^^^*^ is the 
number of electrons from semileptonic decays produced in the signal process that satisfy 
the kinematic requirements using the Monte Carlo information at the generator level. 
At hadron level, the fcy algorithm was applied to all final-state particles with a lifetime 
of r > 0.01 ns and the electron was associated to its parent jet using the generator 
information. 
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All cross sections were measured separately for the two centre-of-mass energies ^/s = 300 
and 318 GeV. Additionally, the cross sections were calculated with the whole data set and 
were corrected to y/s = 318 GeV. The correction factor of ~ 2% was determined with LO 
as well as NLO calculations. 

The charm cross sections were measured using the same procedure. 

9 Systematic uncertainties 

The systematic uncertainties were calculated by varying the analysis procedure and then 
redoing the fit to the likelihood distributions. The following sources were the main con- 
tributors to the systematic uncertainty (the first value in parentheses is the uncertainty 
for beauty, while the second is that for charm): 

• the systematic uncertainty on the description of the dE/dx information was estimated 
by looking at the differences between the various calibration and validation samples. 
Variations in the mean, width and shape of the distributions were evaluated and used 
as a measure of the uncertainty [35]. The resulting uncertainty was found to be (^^% 

/ nyo); 

• the changes in the correction to the distribution in various kinematic ranges were 
taken as a measure of its uncertainty. For = 1.5 GeV the variation was 20% of the 
correction. The changes led to a systematic uncertainty of {'^1% / ^5*^%)- 

In addition, the correction to the charm distribution was varied from zero to that of 
the background sample. This led to an uncertainty of {tl% / 

• a shift of the GAL energy scale in the Monte Garlo simulation by ±3% (±2% / ±5%); 

• reweighting of the direct and non-direct contributions in the Monte Garlo to provide 
a better description of the data (+1% / +3%); 

• the estimated residual number of electrons left in the sample from photon conversions 
as well as from Dalitz decays were varied by 25% and 20% respectively [50]. This led 
to systematic uncertainties of (±1% / =f4%) due to photon conversions and (±1% / 
=F1%) due to Dalitz decays. 

These systematic uncertainties were added in quadrature separately for the negative and 
the positive variations to determine the overall systematic uncertainty of ^9% for the 
beauty and 1^9^% for the charm cross sections. Since no significant dependence of the 
systematic uncertainties on the kinematic variables was observed, the same uncertainty 
was applied to each data point. A 2% overall normalisation uncertainty associated with 
the luminosity measurement was included. 
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A series of further checks were made. The cut on the transverse momentum of the 
electron candidate was varied by ±3%, which is the momentum uncertainty for a track 
with pt = 0.9 GeV. The A0 correction was varied within its uncertainty. The cut on 
AR to associate the decay electron with a jet was varied between 1.5 and 1.0. The 
effect of the gaps between FCAL and BCAL as well as between RCAL and BCAL was 
investigated by varying the cut on the EFO position. Various tests of the signal-extraction 
method were made: e.g. using the likelihood without the Eemc/Eqal or -EcAL/Ptrack 
variables; applying the fit on a signal-enriched sample by making tighter cuts on the input 
variables and varying the fit range. The prior probabilities were recalculated after the 
fit and used as the input for a second fit iteration. Separate fits were made for electron 
and positron candidates for each of the lepton beam particles (e~ and e"*") separately 
as well as for the combined sample. All variations were found to be consistent with 
the expected fluctuations due to statistics; therefore they have not been added to the 
systematic uncertainty. 



10 Theoretical predictions and uncertainties 

QCD predictions at NLO, based on the FMNR programme [51], are compared to the 
data. The programme separately generates processes containing point-like and hadron- 
like photon contributions, which have to be combined to obtain the total cross section. 
The bb and the cc production cross sections were calculated separately. The parton 
distribution functions were taken from CTEQ5M [52] for the proton and GRV-G HO 
[40] for the photon. The heavy-quark masses (pole masses) were set to = 4.75 GeV 
and rric = 1.6 GeV. The strong coupling constant, Aq^j-,, was set to 0.226 GeV. The 
renormalisation, /i/?, and factorisation, fip, scales were chosen to be equal and set to 
fJ'R = fJ'F = sJpt + ^tic) ' where pT is the average transverse momentum of the heavy 
quarks. 

The Peterson fragmentation function [53], with ej, = 0.0035 and ec = 0.035 [54], was used 
to produce beauty and charm hadrons from the heavy quarks. For the bb and cc cross 
sections, the decays into electrons were simulated using decay spectra from Pythia. 

For beauty, both the contributions from prompt and from cascade decays, excluding 
6 — i> r — s> e and b ^ J/ip ^ e~^e~, are taken into account in the effective branching 
fraction. The values were set to 0.221 for the bb and to 0.096 for the cc cross sections [55]. 

For the systematic uncertainty on the theoretical prediction, the masses and scales were 
varied simultaneously to maximise the change in the cross section using the values: 
rri}, = 4.5, 5.0 GeV, rric = 1.35, 1.85 GeV and = ^ip = |\/Pr + "^wc)' ^\/^t + "^Wc)- 
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The effects of different parton density functions as well as variations of e^, within the un- 
certainty of 0.0015 had a small effect on the cross-section predictions and were neglected. 
The parameter ec was varied between 0.02 and 0.07 and the contribution was added in 
quadrature to the systematic uncertainty. The uncertainty on the electron decay spectra, 
evaluated from comparisons to experimental measurements [56,57] and to a simple free- 
quark decay model, was found to be small compared to the total theoretical uncertainty 
and was neglected. 

The uncertainty on the NLO QCD predictions for the total cross section are +25% and 
— 15% for beauty and +45% and —28% for charm. 

The NLO QCD predictions for parton-level jets, reconstructed by applying the fcy al- 
gorithm to the outgoing partons, were corrected for hadronisation effects. A bin-by-bin 
procedure was used according to da = da^i^o " C^^'^, where da^i^o is the cross section for 
partons in the final state of the NLO calculation. The hadronisation correction factor, 
^had^ was defined as the ratio of the dijet cross sections, extracted from the Pythia 
Monte Carlo, after and before the hadronisation process, C^^'^ = da^{Q^"^^ / da^o""^^^ . The 
hadronic corrections were generally small and are given in Tables 1-5. No uncertainty 
was assigned to the correction. 

11 Results 

The visible ep cross sections (quoted at hadron level) for 6-quark and c-quark production 
and the subsequent semileptonic decay to an electron with > 0.9 GeV in the range 
|?7^| < 1.5 in photoproduction events with < 1 GeV^ and 0.2 < y < 0.8 and at least 
two jets with Et > 7(6) GeV, |?7| < 2.5 were determined separately for = 300 GeV 
and y/s = 318 GeV. The measurements are given in Table 1 and are shown in Fig. 9. The 
cross sections at the two different centre-of-mass energies are consistent with each other; 
combining the results leads to a reduced statistical uncertainty. For the complete data 
set (96-00) the cross sections evaluated at i/i = 318 GeV are 

ar = (125 ± ll(stat.)l|?(syst.)) pb, 
af^= (278±33(stat.)l^^(syst.)) pb. 

The visible 6-quark cross section was also determined in the kinematic region of a previous 
ZEUS analysis using muons [5] and is in good agreement. 

The NLO QCD predictions of FMNR (see Table 1) are compared to the data in Fig. 9. 
Good agreement is observed. Also shown in Fig. 9 are expectations of the Pythia Monte 
Carlo. The combined b{b) cross section is a factor 1.75 higher while the c(c) cross section 



11 



is a factor of 1.28 higher than the Pythia prediction (see Section 7). These factors are 
used to scale the Pythia predictions in the following figures. 

Differential cross sections as a function of and ?7e, x°^^, E^r^^ ^ and 77-''^* ^ are shown in 
Figs. 10, 11 and 12, respectively. The variable x'^^ is defined as 

E-pz 

where the sum is over the two highest-energy jets, and corresponds at leading order to the 
fraction of the exchanged-photon momentum in the hard scattering process. The figures 
also show the NLO QCD and the scaled Pythia predictions. The cross-section values 
are given in Tables 2-4. Both the predictions from the NLO QCD calculations as well as 
the scaled Pythia cross sections describe the data well. 

The differential cross sections as a function of the transverse energy of the jet associated 
with the electron from the semileptonic decay, E^^'^^ , were also determined. These cross 
sections are shown in Fig. 13 and given in Table 5. The good agreement with the NLO 
QCD prediction allows the cross section as a function of to be extracted [6]. The 
resulting cross section is shown in Fig. 14 and is also compared with previous ZEUS 
measurements [3,5,6]. The results presented here overlap in with these previous 
measurements and have comparable or smaller uncertainties, giving a consistent picture 
of fe-quark production in ep collisions in the photoproduction regime. 



12 Conclusions 

Beauty and charm production have been measured in dijet photoproduction using semilep- 
tonic decays into electrons. The results were obtained by simultaneously extracting the 
h- and c-quark cross sections using a likelihood ratio optimised for 6-quark production. 
One of the most important variables in the likelihood was the dE/dx measurement in the 
central tracking detector. 

The results were compared to both NLO QCD calculations as well as predictions from 
Monte Carlo models. The NLO QCD predictions are consistent with the data. The Monte 
Carlo models describe well the shape of the differential distributions in the data. The 
results on 6-quark production are also in agreement with a previous less precise ZEUS 
measurement using semileptonic decays into electrons. Within the momentum range 
covered by previous ZEUS measurements using decays into muons, good agreement is 
found. 

The cross section as a function of the transverse momentum of the h quarks has been 
measured over a wide range. The measurements agree well with the previous values. 
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giving a consistent picture of 6-quark production in ep collisions in the photoproduction 
regime, and are well reproduced by the NLO QCD calculations. 
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Table 1: Total cross sections for electrons from h or c quarks in photoproduction 
events, < 1 GeV"^ and 0.2 < y < 0.8, with at least two jets with * > 7(6) GeV , 
\rf^*'\ < 2.5 and the subsequent semileptonic decay to an electron with pfp > 0.9 GeV 
and |r/^| < 1.5. The values are given separately for ^/s = 300 GeV (96 — 97) and 
y/s = 318 GeV (98 — 00) as well as for the complete data set (96 — 00) extrapolated 
to a/s = 318 GeV . The first error is statistical and the second is systematic. In 
addition, the NLO QGD prediction and its uncertainty is given, after applying the 
appropriate hadronisation correction (C^^'^, C^^'^). 
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Table 2: Differential electron cross sections as a function of p^ and rj^ for the 
complete data set. For further details see the caption of Table 1. 
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Table 3: Differential cross sections as a function of x° for the complete data 
set. For further details see the caption of Table 1. 
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12.0±1.9t?;? 




0.97 


28.0±4.7t|^ 


op: + 14 


1.07 


13 : 16 


8.3±1.6l°:? 




1.18 


5.9±3.4+J;° 


11 7+'^-° 

' -2.9 


1.03 


16 : 30 


i.oo±o.38l|^:|]^ 




1.22 


1.5±l.l^°i 


1 + 1.2 


0.89 



^jet 1 


rfcTfe/c/r/j^* ^ 




^had 


dajdrj^"^ ^ 




^had 




(pb) 


(pb) 




(pb) 


(pb) 




-1.0 : -0.25 


24.9±5.2+^:^ 


17.5t|i 


0.82 


73±14+^3 


qq+64 

^^-26 


0.95 


-0.25 : 0.5 


47.6±8.2+|^ 


42.6ti2^7 


1.01 


177±24;?i 




1.05 


0.5 : 1.5 


49.3±7.8l|^ 


30.4;^-9 


0.91 


71±17l^2 


i06ii 


1.04 


1.5 : 2.5 


23.7±5.5t^:? 


q 0+3.6 


0.76 


8±15l| 




1.01 



Table 4: Differential cross sections for the most energetic jet as a function of E; 
and f]^'^^ for the complete data set. For further details see the caption of Table 1. 



Ep"' 


dab/E^-^^^ 




^had 


dac/ E^-^^^ 




^had 


(GeV) 


(pb/ GeV) 


(pb/ GeV) 




(pb/ GeV) 


(pb/ GeV) 




6 : 10 


16.1±1.8t};^ 


i2.3li:^ 


0.67 


42.2±5.2t3:7 


64i38 


1.00 


10 : 15 


6.6±1.3;j|| 




1.00 


22.3±4.2+||^ 


19.6^5^1 


1.06 


15 : 30 


2.1±0.6l°i 


1.08l°:« 


1.21 


0.3±1.9t!]:} 


1 7+1.2 
' -0.5 


0.87 



Table 5: Differential cross sections of E^-^'^ for the jet associated to the electron 
from beauty or charm decays for the complete data set. For further details see the 
caption of Table 1. 
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Figure 1: The mean dE/dx measured in the CTD, (dE/dx), as a function of 
/37 for different samples of identified particles as denoted in the figure. The curve 
shows a physically motivated parametrisation of the {dE/dx) dependence on P'j. 
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Figure 2: Distribution of the relative difference between the observed fciE'/dxobsj 
and predicted (dE / dxpj-ed) specific energy loss for the track sample with ritrunc = 23. 
The inset shows the same distribution with a logarithmic ordinate scale. 



20 



0.5 



1.5 



2 2.5 

dE/dx (mip) 



10" 



0.5 



1.5 



2.5 3 3.5 
,rel , 



Pf (GeV) 




Figure 3: Normalised distributions of the five input variables used in the likelihood 
for the electron candidates extracted from the Monte Carlo, without applying the 
cuts on dE/dx and Eemc/Ecal- The solid line shows the distribution for electrons 
from semileptonic b-quark decays, the dashed line for c-quark decays and the dotted 
line the background (Bkg). 
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Figure 4: The distribution of the likelihood ratio for electron candidates, Nc^nd, 
in data compared to the Monte Carlo expectation after the fit described in the text. 
The arrow indicates the region included in the fit (—2 InT < 10). The shaded areas 
show the fitted contributions from b quarks, c quarks and background as denoted in 
the figure. 
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Figure 5: Distributions of the five input variables of the likelihood for the electron 
candidates used in the fit (—2 InT < 10). All cuts have been applied except dE/dx > 
1.1 in a) and E'emcZ-^'CAL > 0.9 in b) (the cuts are indicated in the figure). The 
shaded areas show the contributions from h quarks, c quarks and background as 
denoted in the figure, after applying the scale factors from the fit. 
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Figure 6: Distributions of the five input variables of the likelihood for a beauty 
enriched sample (candidates with — 21nT < 1). In these plots the ordinate is shown 
on a linear scale. Other details as in the caption of Fig. 5. 
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Figure 7: Distributions of the five input variables of the likelihood for a charm 
enriched sample (candidates with —2 In Te^c < 1-5^). In these plots the ordinate is 
shown on a linear scale. Other details as in the caption of Fig. 5. 
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Figure 8: Distributions of E^-''^*' and rj^''^^ of the jet associated with the electron, 
and of the electron candidate. Figures a)-c) contain all electron candidates 
satisfying —2\nTe,b < 10; d)-f) and g)-i) show the same distributions for the beauty 
(—2lnTg^h < (ii^d charm (—2\nTe^c < l-5j enriched samples, respectively. Other 
details as in the caption of Fig. 5. 
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Figure 9: Total cross sections for electrons from b and c quarks in photoproduction 
events, < 1 GeV'^ and 0.2 < y < 0.8, with at least two jets with Et > 7(6) GeV, 
\rj\ < 2.5 and the subsequent semileptonic decay to an electron with px > 0.9 GeV 
and \rj\ < 1.5. The measurements are shown as points. The inner error bar shows 
the statistical uncertainty and the outer error bar shows the statistical and sys- 
tematic uncertainties added in quadrature. The solid line shows the NLO QGD 
prediction after hadronisation corrections, with the theoretical uncertainties indi- 
cated by the band; the dashed line shows the prediction from Pythia. 
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Figure 10: Differential cross sections as a function of a), c) the transverse 
momentum and h), d) the pseudorapidity of the electrons. Plots a) and b) are for b- 
quark production while c) and d) are for c- quark production. The measurements are 
shown as points. The inner error bar shows the statistical uncertainty and the outer 
error bar shows the statistical and systematic uncertainties added in quadrature. 
The solid line shows the NLO QCD prediction after hadronisation corrections, with 
the theoretical uncertainties indicated by the band; the dashed line shows the scaled 
prediction from Pythia. 
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Figure 11: Differential cross sections as a function of x°^^. a) shows the dis- 
tribution for electrons from b-quark production while b) shows c-quark production. 
Other details as in the caption of Fig. 10. 
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Figure 12: Differential cross sections as a function of a), c) the transverse 
energy and b), d) the pseudorapidity of the highest-energy jet. Plots a) & b) show 
the distributions for electrons from b- quark production while plots c) & d) show 
those for c-quark production. Other details as in the caption of Fig. 10. 
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Figure 13: Differential cross sections for a) b-quark and b) c-quark production 
as a function of the transverse energy of the jet associated to the electron. Other 
details as in the caption of Fig. 10. 
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Figure 14: Differential cross section for h-quark production as a function of 
transverse momentum, p^, compared to the results of previous ZEUS measurements 
as indicated in the figure. The measurements are shown as points. The inner error 
bar shows the statistical uncertainty and the outer error bar shows the statistical 
and systematic uncertainties added in quadrature. The solid line shows the NLO 
QCD prediction from the FMNR program with the theoretical uncertainty shown as 
the shaded band. 
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